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Abstract 

The  matrix  protein  VP40  is  an  indispensable  component  of  viral  assembly  and  budding  by  the  Ebola  virus.  VP40  is  a  monomer  in 
solution,  but  can  fold  into  hexameric  and  octameric  states,  two  oligomeric  conformations  that  play  central  roles  in  the  Ebola  viral 
life  cycle.  While  the  X-ray  structures  of  monomeric  and  octameric  VP40  have  been  determined,  the  structure  of  hexameric  VP40  has 
only  been  solved  by  three-dimensional  electron  microscopy  (EM)  to  a  resolution  of  ~30  A.  In  this  paper,  we  present  the  refinement 
of  the  EM  reconstruction  of  truncated  hexameric  VP40  to  ~20  A  and  the  construction  of  an  all-atom  model  (residues  44-212)  using 
the  EM  model  at  ~20  A  and  the  X-ray  structure  of  monomeric  VP40  as  templates.  The  hexamer  model  suggests  that  the  monomer- 
hexamer  transition  involves  a  conformational  change  in  the  N-terminal  domain  that  is  not  evident  during  octamerization  and  there¬ 
fore,  may  provide  the  basis  for  elucidating  the  biological  function  of  VP40. 
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1.  Introduction 

The  Ebola  virus  is  a  non-segmented,  negative  strand 
RNA  virus  that  causes  hemorrhagic  fever  in  humans 
with  a  mortality  rate  exceeding  70%  (Feldmann  et  al., 
2003).  The  matrix  protein  VP40  is  the  most  abundant 
protein  in  the  Ebola  virus,  and  plays  significant  struc¬ 
tural  and  functional  roles  in  the  Ebola  life  cycle.  Posi- 
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tioned  beneath  the  viral  membrane,  VP40  may  provide 
structural  integrity  to  the  viral  particle  by  linking  the 
membrane  to  the  nucleocapsid  (Feldmann  and  Klenk, 
1996).  In  addition,  VP40  plays  a  central  role  in  viral 
assembly  and  budding  (Aman  et  al.,  2003;  Bavari  et  al., 
2002;  Panchal  et  al.,  2003;  Timmins  et  al.,  2004).  VP40 
binds  and  recruits  specific  human  cellular  proteins 
(Harty  et  al.,  2000;  Licata  et  al.,  2003;  Timmins  et  al., 
2003b)  to  the  budding  site,  allowing  the  virus  to  exploit 
cellular  pathways  in  its  assembly  and  vesicular  release 
(Jasenosky  et  al.,  2001;  Timmins  et  al.,  2001). 

The  X-ray  structure  of  monomeric  VP40  (PDB  entry 
1ES6,  residues  44-321)  has  been  determined  to  a 
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resolution  of  2.0  A  and  reveals  that  VP40  consists  of 
structurally  homologous  N-  and  C-terminal  domains 
(Dessen  et  al.,  2000).  The  two  domains  are  folded  into 
similar  P-sandwiches  composed  of  two  P-sheets  of 
three  antiparallel  strands  with  three  a-helices  packed 
laterally  to  each  P-sandwich  (Fig.  lA).  The  two 
domains  are  oriented  at  an  angle  of  60°  to  one  another, 
are  linked  by  a  large  flexible  loop,  and  are  loosely  asso¬ 
ciated.  In  contrast  to  their  structural  similarity,  the  N- 
and  C-terminal  domains  have  distinct  biological  roles. 
The  N-terminal  domain  is  central  to  oligomerization 
and  the  C-terminal  domain  to  membrane-binding 
(Panchal  et  al.,  2003;  Ruigrok  et  al.,  2000;  Scianima- 
nico  et  al.,  2000). 

While  it  exists  as  a  monomer  in  solution,  VP40  can  also 
fold  into  hexameric  and  octameric  states,  when  treated 
with  urea,  bound  to  liposomes  or  truncated  at  its  C-termi- 


nus  (Ruigrok  et  al.,  2000;  Scianimanico  et  al.,  2000;  Tim¬ 
mins  et  al.,  2003a).  The  X-ray  structure  of  octameric  VP40 
has  been  determined  to  a  resolution  of  1.6  A  and  shows  a 
pore-like  structure  composed  of  four  antiparallel  dimers 
complexed  with  short  RNA  strands  (triribonucleotide  5'- 
U-G-A-3')  (PDB  entry  1H2D,  residues  69-191(192)) 
(Gomis-Ruth  et  al.,  2003)  (Fig.  IB).  With  the  eight  RNA 
segments  bound  at  the  dimer-dimer  interfaces,  the  pore¬ 
like  structure  of  octameric  VP40  has  an  inside  diameter  of 
17  A  (~30  A  in  the  absence  of  the  RNA  strands),  an  out¬ 
side  diameter  of  84  A,  and  a  height  of  42  A  (Fig.  1C). 
Based  upon  analyses  of  the  monomeric  and  octameric 
VP40  X-ray  structures,  it  is  apparent  that  the  monomer- 
octamer  transition  involves:  (1)  the  unfolding  of  the  N- 
terminal  segment  and  (2)  a  conformational  switch  of  the 
N-  and  C-terminal  domains,  in  which  the  C-  terminal 
domain  moves  to  open  up  the  oligomerization  interface. 


Fig.  1.  X-ray  structures  of  monomeric  and  octameric  VP40.  (A)  In  the  X-ray  structure  of  monomeric  VP40  (31-212),  the  N-  and  C-terminal  domains 
are  colored  yellow  and  green,  respectively.  The  N-  and  C-terminal  loops  are  labeled  and  colored  orange  and  purple,  respectively,  to  highlight  their 
relative  positions.  The  N-  and  C-terminal  domains  are  unlinked  because  residues  195-200  were  untraceable  in  the  X-ray  structure.  Furthermore,  the 
N-terminal  segment  (residues  31-43)  was  not  located  on  the  electron  density  map,  presumably  because  of  its  conformational  flexibility.  The  dimen¬ 
sions  of  the  structure  are  40  x  50  x  25  A.  (B  and  C)  The  X-ray  structure  of  octameric  VP40  (69-194)  is  complexed  with  eight  short  RNA  single 
strands  (triribonucleotide  5'-U-G-A-3').  The  eight  VP40  subunits  are  rendered  in  ribbon  and  are  colored  coded.  The  eight  RNA  strands  are  drawn 
as  CPK  models  with  the  carbon  atoms  colored  green,  nitrogen  atoms  blue,  oxygen  atoms  red,  and  phosphorous  atoms  purple.  (B)  The  N-  and  C-ter- 
mini  for  one  antiparallel  dimer  (orange  and  cyan  ribbon)  are  labeled  as  well  as  one  adjacent  subunit  (red  ribbon).  The  pore-like  structure  of 
octameric  VP40  has  an  outside  diameter  of  84  A  and  an  inside  diameter  of  1 7  A  in  the  presence  of  RNA  and  of  30  A  in  the  absence  of  RNA.  (C)  The 
N-  and  C-termini  for  the  subunits  forming  a  dimer-dimer  interface  (cyan  and  red  ribbon)  are  labeled.  A  side  view  shows  that  the  pore-like  structure 
of  octameric  VP40  has  a  height  of  42  A.  For  clarity,  only  four  of  the  eight  VP40  subunits  are  shown. 
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Fig.  2.  EM  reconstruction  of  truncated  (31-212)  and  full-length  (31-326)  VP40.  The  left  column  shows  that  top  view  of  each  EM  reconstruction,  the 
middle  column  was  obtained  by  a  90°  rotation  of  the  top  view  along  the  x-axis,  and  the  right  column  is  a  60°  rotation  of  the  middle  column  along  the 
y-axis.  (A)  The  EM  reconstruction  of  truncated  hexameric  VP40  at  a  resolution  of  ~30  A  reveals  a  pore-like  structure  formed  by  the  oligomerization 
of  the  N-terminal  domain.  The  pore-like  structure  has  an  outside  diameter  of  ~90  A,  an  inside  diameter  of  ~30  A,  and  a  height  of  ~30  A.  The  homo¬ 
dimers  that  comprise  the  pore-like  structure  have  a  width  of  55  A  (adapted  from  Ruigrok  et  al.,  2000).  (B)  The  EM  reconstruction  of  the  nearly  full- 
length  hexameric  VP40  refined  to  ~30  A  shows  that  the  C-terminal  domain  alternates  between  positions  above  and  below  the  pore-like  structure. 
The  C-terminal  domain  is  responsible  for  membrane-binding.  The  lower  electron  density  of  the  C-terminal  domains  relative  to  the  N-terminal 
domains  suggests  that  the  C-terminal  domains  are  conformationally  flexible  (adapted  from  Scianimanico  et  al,  2000).  (C)  The  EM  reconstruction  of 
truncated  hexameric  VP40  refined  to  ~20  A  shows  the  surface  features  of  the  pore-like  structure  in  greater  detail.  The  donut-like  shape  of  hexameric 
VP40  at  ~30  A  decomposes  to  two  triangles  offset  by  ~30°.  Additionally,  the  circular  channel  at  ~30  A  becomes  a  cloverleaf  shape  at  ~20  A.  The 
side  view  of  the  pore-like  structure  shows  a  chair-like  conformation. 

Recent  mutagenesis  data  indicate  that  octameric  VP40  is 
essential  for  the  Ebola  virus  life  cycle  (Hoenen  et  al., 

2005). 

At  present,  the  atomic  structure  of  hexameric  VP40 
has  not  been  determined.  However,  previous  structural 
studies  of  hexameric  VP40  using  electron  microscopy 
have  resulted  in  the  three-dimensional  reconstruction  of 
truncated  (residues  31-212)  and  nearly  full-length  (resi¬ 
dues  31-326)  hexameric  VP40  refined  to  ~30  A  (Figs.  2A 
and  B)  (Ruigrok  et  al.,  2000;  Scianimanico  et  al.,  2000). 

Similar  to  octameric  VP40,  hexameric  VP40  has  a  pore¬ 
like  structure.  From  these  reconstructions,  it  is  apparent 
that  hexameric  VP40  is  a  trimer  of  antiparallel  dimers,  in 
which  the  N-terminal  domains  form  a  pore-like  struc¬ 
ture  and  the  C-terminal  domains  alternate  between  posi¬ 
tions  above  and  below  the  pore-like  structure.  The  lower 
electron  density  of  the  C-terminal  domains  relative  to 
the  N-terminal  domains  suggests  that  the  C-terminal 
domains  are  conformationally  flexible.  One  would 
expect  that  since  hexameric  VP40  consists  of  fewer  N- 
terminal  domains  than  octameric  VP40,  its  pore-like 
structure  would  have  a  smaller  outside  diameter  than 
that  of  octameric  VP40,  assuming  a  similar  conforma¬ 
tion  and  orientation  of  the  N-terminal  domains  in  the 


two  oligomeric  states.  However,  the  pore-like  structure 
of  hexameric  VP40  has  an  outside  diameter  of  ~90A 
and  an  inside  diameter  of  30-35  A,  which  are  the 
approximate  diameters  of  octameric  VP40.  Further¬ 
more,  the  pore-like  structure  of  hexameric  VP40  has  a 
height  of  ~30A,  with  each  antiparallel  dimer  having  a 
width  of  "^55  A.  The  similar  outside  diameters  of  hexa¬ 
meric  and  octameric  VP40  indicate  that  there  exist  major 
structural  differences  between  the  N-terminal  domains 
of  the  two  oligomeric  states.  Unfortunately,  because  the 
atomic  structure  of  hexameric  VP40  has  not  been  deter¬ 
mined,  these  differences  cannot  be  delineated. 

As  hexameric  VP40  appears  to  be  crucial  to  viral 
assembly  and  budding,  determining  its  molecular  struc¬ 
ture  is  of  great  interest.  Here,  we  present  two  models  of 
hexameric  VP40,  a  three-dimensional  EM  reconstruc¬ 
tion  of  truncated  hexameric  VP40  at  ~20A  (31-212) 
which  was  used  in  the  model  building  and  the  all-atom 
model  of  hexameric  VP40  (44-212).  The  model  was  built 
using  the  N-terminal  domain  from  the  monomeric  VP40 
crystal  structure  as  the  initial  structure  and  was  fitted  to 
the  EM  map  of  truncated  hexameric  VP40  at  "^20  A.  The 
two  structures  and  the  conformational  changes  associ¬ 
ated  with  hexamerization  are  discussed  herein. 
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2.  Materials  and  methods 

2.7.  General  methods 

A  Silicon  Graphics  Octane  2  workstation  with  Insight 
II  (Accelrys,  San  Diego,  CA)  was  used  to  build  and  visu¬ 
alize  the  models.  All  simulations  were  performed  using 
the  Discover  3.0  program  (Accelrys,  San  Diego,  CA) 
with  the  molecular  mechanics  potentials  set  to  the  amber 
force  field  and  the  non-bonded  interaction  limited  to 
1 3  A.  The  tethered  minimization  and  constrained  molec¬ 
ular  dynamics  protocols  used  to  construct  the  hexameric 
VP40  model  were  previously  described  (Nguyen,  2004). 
The  quality  of  the  models  was  assessed  using  Procheck 
(Laskowski  et  al.,  1993)  and  HINT  (eduSoft,  Richmond, 
VA).  The  root  mean  square  deviation  (r.m.s.  deviation) 
values  were  calculated  relative  to  the  mean  structure 
using  MOLMOL  (Koradi  et  al.,  1996),  and  the  figures 
were  rendered  using  Molscript  (Kraulis,  1991),  Pymol 
(DeLano  Scientific,  San  Carlos,  CA),  and  Insight  II. 

2.2.  EM  model  and  fit 

The  EM  model  was  recalculated  from  previous  data 
with  the  same  projection  matching  method  as  in  Rui- 
grok  et  al.  (2000)  but  by  correcting  the  images  for  the 
contrast  transfer  function  (CTF)  according  to  the  Con¬ 
way  and  Steven  (1999)  method.  The  resolution  of  the 
final  map  is  ~20  A,  as  tested  by  Fourier  shell  correlation. 
The  fit  between  the  EM  map  and  the  molecular  model 
was  performed  manually  using  O,  and  the  result  was 
then  submitted  to  Situs  (Wriggers  and  Birmanns,  2001), 
which  gave  a  correlation  of  0.208. 

2.3.  All-atom  model 

Fig.  3  outlines  the  protocol  used  to  generate  the  all¬ 
atom  model.  The  initial  structure  is  the  monomeric  VP40 
crystal  structure,  which  consists  of  the  N-  and  C-termi- 
nal  domains.  Residues  213-321  of  the  C-terminal 
domain  were  removed.  The  result  is  a  truncated  N-termi- 
nal  domain  (residues  44-212).  Using  the  Ca  trace  of  the 
P-sandwiches  as  guides,  two  copies  of  the  truncated 
monomeric  form  of  the  N-terminal  domain  were  manu¬ 
ally  superimposed  onto  the  antiparallel  dimer  of  the 
octameric  VP40  crystal  structure.  The  protomer-pro- 
tomer  interface  was  both  energetically  and  hydropathi- 
cally  refined.  This  involved:  (1)  manual  adjustment  of 
the  torsion  angles  of  the  side  chains  at  the  protomer- 
protomer  interface  to  relieve  van  der  Waals  violations  of 
greater  than  0.25  A,  (2)  energy  minimization  involving 
up  to  5000  steps  of  Fletcher-Powell  optimization  until 
the  norm  of  the  gradient  was  <  1 .0  kcal/mol  A^.  During 
this  energy  minimization,  the  backbone  atoms  of  the  two 
protomers  were  fixed  in  Cartesian  space  so  that  the  over¬ 
all  structure  of  the  complex  was  maintained,  (3)  hydro- 


Fig.  3.  Schematic  detailing  the  methodology  used  to  model  the  hexa¬ 
meric  VP40  structure.  The  initial  structure  is  the  X-ray  structure  of 
monomeric  VP40  (44—321),  which  consists  of  the  two  structurally 
homologous  N-  and  C-terminal  domains,  (a)  Residues  of  the  C-terminal 
domain  (213-321)  are  removed  and  residues  which  were  missing  in  the 
N-terminal  domain  were  modeled  to  produce  a  N-terminal  domain  with 
a  sequence  of  44-212.  (b)  Using  the  antiparallel  dimers  of  the  octameric 
VP40  crystal  structure  as  a  template,  two  copies  of  the  modeled  N-ter¬ 
minal  domain  were  dimerized  in  an  antiparallel  fashion.  The  protein- 
protein  interface  was  refined  to  the  hydropathic  quality  evident  in  the 
octameric  VP40  crystal  structure,  (c)  Using  the  volume  and  dimensions 
of  the  EM  reconstruction  of  hexameric  VP40  as  restraints,  three  copies 
of  the  model  VP40  dimer  were  positioned  into  a  hexameric  assembly,  (d) 
Because  there  are  large  gaps  in  the  dimer-dimer  interfaces,  it  was  neces¬ 
sary  to  use  constrained  molecular  dynamics  to  produce  tightly  packed 
protein-protein  surfaces.  This  mostly  involved  the  rigid-body  shift  of 
the  |3-sheet  consisting  of  strands  |33,  |34,  and  |35  towards  the  center  of  the 
protein-protein  gap.  One  dimer-dimer  interface  is  refined,  (e  and  f) 
Instead  of  modeling  the  conformational  change  at  the  two  remaining 
dimer-dimer  interfaces,  the  two  subunits  of  the  refined  dimer-dimer 
interface  are  hexamerized.  The  protomer-protomer  interfaces  are  ener¬ 
getically  and  hydropathically  refined,  (g)  The  C-terminal  loop  residues 
which  were  in  the  nearly  extended  conformation  in  the  monomeric 
VP40  crystal  structure  are  packed  using  energy  minimization  and  con¬ 
strained  molecular  dynamics  to  produce  the  final  model. 

pathic  analysis  of  the  refined  protomer-protomer 
interface  using  the  HINT  program  and  comparison  to 
the  hydropathic  quality  of  the  octameric  VP40  crystal 
structure,  and  (4)  constrained  energy  minimization  in 
which  unfavorable  hydrophobic-polar  and  basic-basic 
interactions  as  defined  by  HINT  values  higher  than 
those  in  the  octameric  VP40  crystal  structure  were 
relieved  by  distance  constraints  which  incrementally  sep¬ 
arated  the  atoms  involved  in  the  unfavorable  interaction. 

Initially,  three  of  the  modeled  antiparallel  dimers 
were  arrayed  into  a  ring-like  structure  with  a  diameter  of 
~90A.  Because  of  their  size  and  shape,  the  antiparallel 
dimers  in  this  initial  ring-like  structure  were  not  tightly 
packed,  but  were  separated  by  ~15  A  gaps,  as  measured 
between  Ca  atoms.  The  protomers  were  rotated  in  a 
rigid-body  fashion  in  an  attempt  to  close  these  large 
protein-protein  gaps.  However,  this  approach  was 
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unsuccessful.  Structural  analyses  indicated  that  closing 
these  large  gaps  required  a  major  conformational  change 
in  the  sheet  structure  formed  by  strands  P3,  P4,  and  P5.  If 
these  strands  from  each  protomer  were  moved  in  a  con¬ 
certed  fashion  towards  the  opening  at  the  center  of  the 
dimer-dimer  interface  and  were  packed  against  one 
another  in  an  energetically  and  hydropathically  reason¬ 
able  manner,  the  result  would  be  a  pore-like  structure 
with  an  outside  diameter  of  ~90  A  and  an  inside  diame¬ 
ter  of  "^30  A.  Absent  this  conformational  change  in  the 
tertiary  structure  of  VP40,  the  large  gaps  at  the  dimer- 
dimer  interfaces  could  be  closed  by  the  rigid-body  trans¬ 
lation  of  the  three  antiparallel  dimers  towards  the  center 
of  the  ring-like  structure.  However,  this  would  produce  a 
ring-like  structure  with  an  outside  diameter  of  only 
"^70  A  and  lacking  a  channel  (data  not  shown). 

A  constrained  molecular  dynamics  protocol  was  used 
to  produce  the  necessary  conformational  changes  in 
strands  P3,  P4,  and  pS.  Since  the  antiparallel  dimers  in 
hexameric  VP40  are  related  by  a  threefold  rotational 
axis,  an  efficient  molecular  design  strategy  and  the  one 
which  was  used  here  is  to  close  the  spatial  gap  at  one 
dimer-dimer  interface  and  then  to  copy  the  closed 
dimer-dimer  subunit  to  the  other  two  positions  to  form 
the  closed  hexamer.  Because  the  hexamer  must  be  reas¬ 
sembled  from  different  dimer-dimer  subunits,  it  was  nec¬ 
essary  to  fix  the  positions  of  some  of  the  atoms  of  the 
dimer-dimer  subunits  during  molecular  dynamics  in 
order  to  properly  reassemble  the  hexamer.  We  chose  to 
fix  the  peptide  backbone  of  strands  pi,  P2,  and  P6,  since 
this  prevents  unnecessary  translation  or  rotation  of  the 
protomers  and  additionally  maintains  the  structural 
integrity  of  the  protomer-protomer  interfaces.  During 
the  simulations,  distance  constraints  were  used  to  main¬ 
tain  the  secondary  structure  of  the  protein  and  allowed 
strands  P3,  P4,  and  P5  to  move  as  a  single  unit.  Two  pro¬ 
tomers  at  a  dimer-dimer  interface  were  selected  and  a 
methane  molecule  was  positioned  near  the  center  of  the 
dimer-dimer  interface  but  outside  of  the  ring-like  struc¬ 
ture.  Distance  constraints  were  created  between  the 
methane  carbon  and  the  Lysl27  Ca  atoms  of  the  two 
protomers  at  the  dimer-dimer  interface.  During  the 
molecular  dynamics  simulations,  the  upper  limits  of 
these  distance  constraints  were  incrementally  decreased 
to  bring  the  two  Lysl27  Ca  atoms  closer  to  the  methane 
molecule,  and  consequently,  the  two  sheet  structures 
towards  the  center  of  the  dimer-dimer  interface.  Once 
the  two  sheet  structures  were  brought  within  van  der 
Waals  contact  distance  to  one  another,  tethered  minimi¬ 
zation,  molecular  dynamics,  and  torsional  angle  adjust¬ 
ments  were  used  to  energetically  and  hydropathically 
optimize  the  protein-protein  interactions. 

Using  the  initial  ring-like  structure  as  a  template, 
three  of  the  modeled  dimer-dimer  subunits  were  arrayed 
into  a  hexamer.  With  the  peptide  backbone  of  strands 
pi,  P2,  and  P6  no  longer  fixed  in  Cartesian  space,  the 


protomer-protomer  interfaces  in  the  above  hexameric 
structure  were  refined  using  tethered  minimization, 
molecular  dynamics,  and  torsional  adjustment.  In  the 
final  stage,  residues  187-212,  which  are  part  of  the  large 
flexible  loop  connecting  the  N-  and  C-terminal  domains, 
were  folded  using  constrained  molecular  dynamics  and 
positioned  near  the  outside  corners  of  the  pore-like 
structure.  The  quality  of  the  model  was  assessed  using 
Procheck,  which  showed  98.3%  of  the  residues  in  the 
most  favored  or  additionally  allowed  regions. 

3.  Results  and  discussion 

3.1.  EM  map  refined  to  ^20  A 

The  three-dimensional  EM  reconstruction  of  truncated 
and  full-length  hexameric  VP40  had  previously  been 
reported  at  ~30A  (Ruigrok  et  al.,  2000;  Scianimanico 
et  al.,  2000).  Both  EM  reconstructions  show  that  the  pore¬ 
like  structure  of  hexameric  VP40  has  a  donut-like  shape 
with  a  circular  channel  and  rounded  edges  (Figs.  2 A  and 
B),  which  from  the  side  view,  has  a  chair-like  conforma¬ 
tion.  Here,  refinement  of  the  EM  map  of  truncated  hexa¬ 
meric  VP40  to  ~20A  reveals  some  important  structural 
characteristics.  First,  the  donut-like  structure  evident  in 
the  top  view  at  ~30A  decomposes  into  two  triangular 
structures  at  ~20  A  (Fig.  2C).  The  two  triangles  are  offset 
from  one  another  by  ~30°  relative  to  the  channel  axis. 
Second,  at  the  higher  resolution,  the  channel  of  hexameric 
VP40  adopts  a  cloverleaf-like  shape  (Fig.  2C).  The  pro¬ 
tomer-protomer  interfaces  outline  the  three  leaves  of  the 
cloverleaf  shape,  while  the  dimer-dimer  interfaces  form 
the  bulges  in  the  channel.  The  side  view  of  the  EM  recon¬ 
struction  shows  a  similar  chair  conformation  at  ~20A 
and  at  ~30  A.  As  evident  in  Fig.  2,  more  of  the  structural 
details  of  hexameric  VP40  are  revealed  at  ~20  A  and  these 
features  were  incorporated  into  the  model  building. 

3.2.  Overall  strueture  of  the  hexamer  model 

As  shown  in  Fig.  4,  the  overall  structure  of  the  hexamer 
model  is  strikingly  similar  to  the  EM  reconstruction  at 
~20A.  Consistent  with  EM  reconstructions,  the  atomic 
model  has  an  outside  diameter  of  90  A,  an  inside  diameter 
of  30-35  A,  and  a  height  of  30^0  A.  The  height  decreases 
to  30  A  near  the  dimer-dimer  interface  with  the  dimers 
having  a  width  of  '--SS  A.  Besides  their  similar  dimensions, 
the  hexamer  model  features  a  channel  with  the  distinctive 
cloverleaf  shape  evident  in  the  EM  maps  at  ~20A.  As 
expected,  while  the  leaves  of  the  cloverleaf  are  shaped  by 
the  protomer-protomer  interface,  the  bulges  in  the  chan¬ 
nel  are  formed  by  the  dimer-dimer  interfaces.  The  surface 
potential  at  the  top  of  the  hexamer  model  is  mostly  acidic 
near  the  protomer-protomer  interfaces  and  basic  near  the 
dimer-dimer  interfaces  (Fig.  4A).  From  the  side  view,  the 
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Fig.  4.  Molecular  model  of  truncated  hexameric  VP40  (44-212)  and  its  fit  to  the  EM  map.  The  hexamer  model  is  rendered  in  ribbon  and  the  proto- 
mers  are  colored  coded.  The  N-  and  C-termini  are  shown  for  one  subunit  (rendered  in  green).  Surface  renderings  are  colored  according  to  the  nega¬ 
tive  (red)  and  positive  (blue)  electrostatic  surface  potential.  (A)  The  top  view  shows  that  the  hexamer  model  has  a  triangular  shape  with  an  outside 
diameter  of  90  A  and  an  inside  diameter  of  ~30  A.  The  dimer-dimer  interfaces  are  mostly  basic  and  the  protomer-protomer  interfaces  are  mostly 
acidic.  The  cloverleaf  shape  of  the  channel  is  also  evident.  (B)  A  90°  rotation  along  the  x-axis  of  the  top  view  shows  that  the  model  has  a  height  of  30- 
40  A.  The  height  varies  and  narrows  to  ~30  A  at  the  dimer-dimer  interfaces.  Near  the  dimer-dimer  interfaces,  there  are  both  positive  and  negative 
charge  surfaces.  (C)  A  60°  rotation  along  the  y-axis  of  (B)  shows  that  the  protomer-protomer  interfaces  are  mostly  acidic.  (D)  The  hexamer  model  is 
fitted  to  the  EM  map  at  ~20  A,  resulting  in  a  coefficient  of  0.208.  The  top  view  shows  that  the  model  has  the  triangular  shape  of  the  EM  reconstruc¬ 
tion  as  well  as  the  cloverleaf-like  channel.  The  side  view  (a  90°  rotation  along  the  x-axis)  shows  the  chair-like  conformation  of  hexameric  VP40. 
While  much  of  the  hexamer  model  fits  into  the  EM  map,  the  C-terminal  residues  190-212  are  outside  of  the  electron  density.  Residues  190-212  form 
the  large  flexible  loop  connecting  the  N-  and  C-terminal  domains  and,  therefore,  are  highly  disordered. 


surface  potential  is  both  positive  and  negative  at  the 
dimer-dimer  interfaces  (Fig.  4B),  but  mostly  negative  at 
the  protomer-protomer  interfaces  (Fig.  4C). 

Each  protomer  in  the  hexamer  model  is  composed  of 
two  P-sheets  consisting  of  three  antiparallel  strands  and 
three  a-helices.  A  similar  secondary  structure  is  evident 
for  the  monomeric  and  the  octameric  states  of  the  N-ter- 
minal  domain.  As  a  consequence  of  the  approach  used  to 
generate  the  model,  the  protomers  of  the  model  are  not 
identical,  but  have  a  r.m.s.  deviation  of  ~lA  for  the 
backbone  atoms  of  residues  44-194. 

3.3.  Docking  of  the  model  into  the  EM  map 

To  validate  the  hexamer  model,  it  was  docked  into  the 
EM  map  of  truncated  hexameric  VP40  at  ~20A.  The 


result  was  a  correlation  of  0.208.  Fig.  4D  shows  that 
while  much  of  the  hexamer  model  fits  into  the  EM  map 
at  ~20  A,  the  C-terminal  segments  (residues  190-212)  of 
the  hexamer  model  are  outside  of  the  electron  density. 
Residues  190-212  form  the  large  loop  connecting  the  N- 
and  C-terminal  domains  in  monomeric  VP40  and  are 
highly  disordered,  as  indicated  by  the  fact  that  residues 
195-200  were  untraceable  in  the  monomeric  VP40  X-ray 
structure.  Because  of  their  conformational  flexibility, 
residues  190-212  may  not  be  detectable  in  the  EM 
images.  Similarly,  the  N-terminus  may  be  highly  disor¬ 
dered  (residues  31-43  were  not  located  in  the  electron 
density  map  of  monomeric  VP40)  and  not  visible  by  EM. 
Although  it  is  clear  that  the  unfolding  of  the  N-terminal 
segment  is  a  critical  step  in  VP40  oligomerization,  there 
is  little  structural  information  to  suggest  the  conforma- 
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Fig.  5.  Stereoview  of  the  superimposition  of  the  calculated  average 
protomer  of  the  hexamer  model  (44-194)  and  the  N-terminal  domain 
of  the  monomeric  VP40  X-ray  structure  (44-194).  The  average  pro¬ 
tomer  of  the  hexamer  model  is  rendered  in  yellow  ribbon  and  the 
monomer  crystal  structure  in  cyan  ribbon.  The  two  structures  were 
superimposed  along  the  Ca  trace  of  strands  pi,  p2,  and  p6.  (A)  While 
much  of  the  two  structures  can  be  overlapped,  there  are  considerable 
differences  between  the  two,  particularly  for  strands  p3,  p4,  and  p5 
which  form  a  three-stranded  p-sheet.  (B)  A  rotation  of  90°  along  the 
x-axis  shows  the  top  view  of  the  superimposition  and  illustrates  the 
different  conformational  space  occupied  by  the  N-terminal  domains  of 
the  hexamer  model  and  the  monomer  crystal  structure. 

tion  and  position  of  the  N-terminal  segments  in  hexa- 
meric  VP40.  Accordingly,  while  the  construct  used  in  the 
EM  studies  consisted  of  residues  31-212,  residues  31^3 


were  not  included  in  the  hexameric  VP40  model.  The  N- 
terminal  segments  of  the  hexameric  VP40  model  were 
positioned  at  the  dimer-dimer  interface  and  on  the  out¬ 
side  of  the  pore-like  structure  (Fig.  4C).  Finally,  the  ring 
interior  of  the  hexamer  model  has  more  of  a  triskelion 
appearance.  This  could  easily  be  due  to  the  use  of  stain 
that  induces  a  Fresnel  fringe  at  the  inside  of  the  particle. 

3.4.  Comparison  to  the  monomeric  VP40  X-ray  structure 

The  N-terminal  domains  of  the  hexamer  model  and 
of  the  monomeric  VP40  crystal  structure  possess  the 
same  secondary  structural  elements,  but  have  distinct 
tertiary  structures.  Figs.  5 A  and  B  show  the  superimposi¬ 
tion  of  the  N-terminal  domain  (residues  44-194)  from 
the  X-ray  structure  of  monomeric  VP40  and  the  calcu¬ 
lated  average  protomer  from  the  hexamer  model.  The 
superimposition  shows  that  while  the  two  N-terminal 
domains  occupy  some  common  conformational  space, 
major  structural  differences  exist  between  the  two. 

Based  on  the  superposition  of  strands  (31,  (32,  and  P6, 
the  two  structures  have  a  backbone  r.m.s.  deviation  of 
2.9  A  for  residues  44—194  (Fig.  6).  This  relatively  low 
r.m.s.  deviation  value  belies  the  major  conformational 
changes  in  the  N-terminal  domain  as  a  result  of  hexa- 
merization.  While  strands  pi,  P2,  and  P6,  which  form 
one  antiparallel  P-sheet  in  the  N-terminal  domain,  have 
an  average  backbone  r.m.s.  deviation  of  only  0.6  A, 
strands  P3,  P4,  and  P5,  which  make  up  the  second 
P-sheet  structure  in  the  N-terminal  domain,  have  an 
average  backbone  r.m.s.  deviation  of  4.6  A.  The  signifi¬ 
cantly  higher  r.m.s.  deviation  of  strands  P3,  P4,  and  P5  is 
not  unexpected,  given  that  these  strands  were  used  to 
close  the  protein-protein  gaps  during  model  building. 
Because  of  the  modeling  methodology  used  to  build  the 
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Fig.  6.  Backbone  r.m.s.  deviation  values  of  the  calculated  average  protomer  of  the  hexamer  model  compared  to  the  N-terminal  domains  in  the  X-ray 
structures  of  monomeric  (44-194)  and  octameric  VP40  (69-191).  The  squares  (■)  depict  the  r.m.s.  deviation  values  between  the  calculated  average 
hexamer  protomer  and  monomeric  form  of  the  N-terminal  domain  and  the  triangles  (A)  between  the  calculated  average  hexamer  protomer  and  the 
octameric  state  of  the  N-terminal  domain.  Because  the  monomeric  and  octameric  states  of  the  N-terminal  domain  are  structurally  similar,  the  two 
sets  of  r.m.s.  deviation  values  parallel  one  another.  The  highest  r.m.s.  deviation  values  are  observed  for  residues  127-129,  which  form  a  turn  connect¬ 
ing  strands  |33  and  |34. 
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hexameric  VP40  model,  the  turn  residues  127-129 
underwent  the  largest  conformational  changes  and  have 
an  average  r.m.s.  deviation  value  of  11.6  A.  The  model 
suggests  that  the  monomer-hexamer  transition  involves 
a  shift  in  two  P-sheet  structures  of  the  N-terminal 
domain.  As  a  measure  of  this  shift,  the  centers  of  the  two 
P-sheets  move  from  a  separation  of  14.3  A  in  the  mono¬ 
meric  VP40  crystal  structure  to  16.6  A  in  the  hexamer 
model,  a  translation  of  2.3  A. 

3.5.  Comparison  to  the  octameric  VP40  X-ray  structure 

Although  hexameric  VP40  consists  of  fewer  subunits 
than  octameric  VP40,  the  pore-like  structures  of  the  hex¬ 
americ  VP40  model  and  the  octameric  VP40  crystal 
structure  have  similar  outside  diameters  of  90  and  84  A, 
respectively,  an  indication  that  the  N-terminal  domains 
that  make  up  the  two  pore-like  structure  have  different 
conformations.  Fig.  7  shows  a  side-by-side  view  of  the 
antiparallel  dimers  from  the  hexamer  model  and  the 
octamer  X-ray  structure.  While  consisting  of  similar  sec¬ 
ondary  structures,  the  tertiary  structures  of  the  two 
dimers  are  strikingly  different.  As  illustrated  in  Fig.  7B, 


Fig.  7.  Structural  comparison  of  the  antiparallel  dimers  in  the  hexa¬ 
meric  VP40  model  (residues  70-190)  and  in  the  octameric  VP40  X-ray 
structure  (70-190).  Both  dimers  are  rendered  in  ribbon.  The  two  sub¬ 
units  of  the  antiparallel  dimer  in  the  hexamer  model  are  colored  blue 
and  cyan  (left),  while  those  of  the  octamer  crystal  structure  are  colored 
green  and  orange  (right).  The  two  dimers  are  positioned  with  a  similar 
orientation  and  illustrated  side-by-side.  The  sheet  structures  in  one 
subunit  of  the  hexamer  model  and  in  one  subunit  of  the  octamer  crys¬ 
tal  structure  are  labeled.  (A)  While  consisting  of  similar  secondary 
structure,  the  tertiary  structures  of  the  two  dimers  are  markedly  differ¬ 
ent.  (B)  A  90°  rotation  of  the  two  structures  along  the  x-axis  shows 
that  the  dimer  from  the  hexamer  model  is  a  more  open  structure. 


the  antiparallel  dimer  of  the  hexamer  model  has  a  more 
open  structure  than  that  of  the  octamer. 

Based  on  the  superimposition  of  strands  pi,  P2,  and 
P6,  subunit  A  of  the  octameric  VP40  crystal  structure 
and  the  calculated  average  protomer  of  the  hexamer 
model  have  a  backbone  r.m.s.  deviation  of  3.5  A  for  resi¬ 
dues  69-190  (Fig.  6).  Additionally,  there  are  large  differ¬ 
ences  in  the  P-sandwich  structures  of  the  hexamer  model 
and  octamer  crystal  structure.  While  the  P-sheet  formed 
by  strands  pi,  P2,  and  P6  has  an  average  r.m.s.  deviation 
of  0.9  A,  the  second  P-sheet  composed  of  strands  P3,  P4, 
and  P5  has  an  average  value  of  5.6  A.  Similar  to  the 
monomer-hexamer  comparison,  the  highest  values  are 
observed  for  residues  127-129  which  form  a  turn  con¬ 
necting  strands  P3  and  P4.  In  the  octameric  VP40  crystal 
structure,  residues  127-129  are  not  hydrogen  bonded 
and  are  positioned  near  the  center  of  the  channel.  In  con¬ 
trast,  the  hexamer  model  suggests  that  the  127-129  loop 
is  positioned  at  the  dimer-dimer  interface  in  hexameric 
VP40.  These  r.m.s.  deviation  values  indicate  that  the  two 
P-sheets  in  VP40  move  farther  apart  during  hexameriza- 
tion  than  during  octamerization.  In  the  octamer  X-ray 
structure  and  the  hexamer  model,  the  centers  of  the  two 
P-sheets  are  separated  by  14.7  and  16.6  A,  respectively. 

3.6.  Comparison  of  the  protein-protein  interfaces 

Because  the  boundaries  between  the  VP40  subunits  in 
the  EM  reconstructions  are  not  well  defined,  the  features 
of  the  protein-protein  interfaces  in  hexameric  VP40  are 
largely  unknown.  By  providing  atomic  detail  to  the  EM 
reconstructions,  the  hexamer  model  suggests  potential 
interactions  at  the  protein-protein  interfaces  and  may 
provide  the  basis  for  future  work. 

There  are  two  types  of  protein-protein  interfaces  in 
the  hexamer  model  and  the  octamer  crystal  structure. 
The  first,  the  protomer-protomer  interface,  reflects  the 
contacts  formed  between  the  protomers  within  the 
dimers  of  the  two  oligomeric  conformations  of  VP40, 
while  the  second,  the  dimer-dimer  interface,  reflects  the 
interactions  between  dimers.  The  hexamer  model  sug¬ 
gests  that  hexameric  and  octameric  VP40  have  similar 
protomer-protomer  interfaces  and  strikingly  different 
dimer-dimer  interfaces. 

In  the  octameric  VP40  crystal  structure,  hydrophobic 
contacts  dominate  the  protomer-protomer  interfaces 
with  Ile74,  Trp95,  Pro97,  Phel61,  and  Ilel82  forming  a 
hydrophobic  core.  Similarly,  in  the  hexamer  model, 
hydrophobic  contacts  comprised  of  Ile74,  Trp95,  Pro97, 
and  He  182  appear  to  stabilize  the  protomer-protomer 
interfaces  (Fig.  8A).  Additionally,  in  octameric  VP40, 
this  hydrophobic  core  is  complemented  by  salt  bridge 
interactions  from  Glul60  to  Argl48  and  Argl51.  The 
hexamer  model  suggests  that  this  Argl48-Glul60  salt 
bridge  is  replaced  by  a  cation-Ti  interaction  (Flocco  and 
Mowbray,  1994;  Mitchell  et  al.,  1994)  between  the  ring 
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Fig.  8.  Notable  interactions  at  the  protein-protein  interfaces  in  the 
hexamer  model.  The  two  interface  subunits  are  rendered  as  ribbon  and 
colored  either  green  or  cyan.  The  identified  residues  are  rendered  in 
sticks  with  the  carbon  atoms  colored  white,  the  nitrogen  atoms  blue, 
and  the  oxygen  atoms  red.  (A)  The  protomer-protomer  interface  is 
dominated  by  hydrophobic  contacts.  Ile74,  Trp95,  Pro97,  and  Ilel82 
form  a  hydrophobic  core.  Complementing  this  hydrophobic  core  is  a 
cation-TC  interaction  between  Trp95  and  Argl48  of  the  neighboring 
subunit.  (B)  The  dimer-dimer  interfaces  also  feature  a  cation-Ti  inter¬ 
action,  but  in  this  instance  between  Phel25  and  Lysl27.  (C)  The  bulges 
evident  in  the  channel  interior  of  the  model  and  the  EM  reconstruc¬ 
tion  of  truncated  hexameric  VP40  at  20  A  occur  at  the  dimer-dimer 
interfaces  and  are  mostly  hydrophobic,  consisting  of  Leu  158,  Phel61, 
Vall62,  Prol64,  Prol65,  and  Vall66. 


of  Trp95  and  the  guanidinium  group  of  Argl48  from 
neighboring  subunits  of  hexameric  VP40. 

The  hexamer  model  suggests  that  hydrophobic  resi¬ 
dues  are  critical  to  the  stability  of  the  dimer-dimer  inter¬ 
faces  in  hexameric  VP40  with  the  aromatic  ring  of 
Phel25  playing  a  central  role.  In  addition  to  providing 
hydrophobic  contacts,  the  ring  of  Phel25  forms  a  cat- 
ion-TT  interaction  with  the  side  chain  amino  group  of 
Lysl27  (Fig.  8B).  In  addition,  clusters  of  hydrophobic 
residues  line  the  channel  interior  at  the  dimer-dimer 
interfaces.  These  clusters  form  the  three  bulges  which  are 
evident  in  the  EM  reconstruction  at  ~20A  and  are 
responsible  for  the  cloverleaf-like  shape  of  the  channel. 
Leul32,  Phel57,  Leul58,  Phel61,  Vall62,  Prol64, 
Prol65,  and  Vall66  comprise  these  clusters  (Fig.  8C). 

Although  it  has  been  shown  that  RNA  stabilizes  hex¬ 
americ  VP40  in  vitro,  in  this  paper,  the  hexameric  VP40 
model  was  constructed  in  the  absence  of  RNA.  The  pres¬ 
ence  of  RNA  would  not  significantly  alter  the  overall 
structure  of  the  hexameric  VP40  model,  but  would 
necessitate  the  conformational  change  of  certain  resi¬ 
dues  in  order  for  hexameric  VP40  to  bind  RNA.  In  the 
X-ray  structure  of  octameric  VP40,  side  chains  of 
Phel25  and  Argl34  of  VP40  form  strong  interactions 
with  the  RNA  strands;  however,  because  of  the  absence 
of  RNA,  side  chains  of  Phel25  and  Argl34  in  the  hexa¬ 
mer  model  form  residue-residue  interactions.  As  shown 
in  Fig.  8B,  the  hexameric  VP40  model  suggests  that  the 
aromatic  ring  of  Phel25  forms  a  cation-Ti  interaction 
with  the  side  chain  amino  group  of  Lysl27.  Addition¬ 
ally,  the  side  chain  of  Argl34  forms  multiple  strong 
hydrogen  bonds.  While  it  may  be  possible  to  model  the 
hexameric  VP40-RNA  complex,  the  unspecific  binding 
of  RNA  to  hexameric  VP40  and  the  conformational 
flexibility  of  the  RNA  backbone  makes  modeling  the 
hexamer-RNA  complex  extremely  difficult.  However,  as 
shown  in  Fig.  4,  the  hexameric  VP40  model  features 
large  patches  of  positive  electrostatic  surface  potential  at 
the  dimer-dimer  interfaces,  where  RNA  could  bind  non- 
specifically  and  provide  additional  stability  to  the  pro¬ 
tein-protein  interfaces  as  suggested  by  biochemical  data 
(Timmins  et  al.,  2003a).  In  contrast,  RNA  binding  in 
octameric  VP40  is  specific  and  plays  a  structural  role  in 
stabilizing  the  protein-protein  interfaces. 

3. 7.  Conformational  states  of  VP40 

Besides  the  monomeric  form,  VP40  can  adopt 
dimeric,  hexameric,  and  octameric  states  (Timmins  et  al., 
2003a).  Assuming  that  dimeric  VP40  is  the  building 
block  for  both  hexameric  and  octameric  VP40,  why  can 
VP40  oligomerize  into  both  higher  order  states?  A  struc¬ 
tural  analysis  of  the  hexamer  model  and  the  octamer 
crystal  structure  reveals  that  the  antiparallel  dimers  that 
comprise  both  structures  mostly  differ  in  the  spatial  sep¬ 
aration  of  their  two  P-sheets.  The  degree  of  twist  that 
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can  be  achieved  between  the  two  P-sheets  of  the  N-ter- 
minal  domain  may  dictate  which  of  the  higher  order 
states  VP40  forms.  The  centers  of  the  two  P-sheets  are 
A  farther  apart  in  the  hexamer  model  than  they  are 
in  the  octamer  crystal  structure.  This  simple  conforma¬ 
tional  switch  in  the  relative  position  of  the  two  P-sheets 
of  the  N-terminal  domain  may  allow  VP40  to  assume 
two  different  oligomeric  states  and  as  a  corollary,  per¬ 
form  multiple  functional  and  structural  roles. 

4.  Conclusion 

Although  the  atomic  structure  of  hexameric  VP40  has 
so  far  eluded  determination,  we  believe  that  there  is  suffi¬ 
cient  biochemical  and  structural  data  on  VP40  to  build  a 
reasonable  molecular  model  of  the  pore-like  structure  of 
hexameric  VP40.  Since  hexameric  VP40  plays  an 
important  role  in  viral  assembly  and  budding  by  the 
Ebola  virus,  its  atomic  structure  may  prove  useful  in  elu¬ 
cidating  its  function.  Here,  the  two  models  are  presented. 
The  first  is  an  EM  reconstruction  of  truncated  hexameric 
VP40  at  "^20  A  and  the  second  is  an  all-atom  model  which 
was  fitted  to  the  EM  map  with  a  correlation  of  0.208.  The 
molecular  model  indicates  that  hexamerization  involves  a 
major  conformational  change  in  the  N-terminal  domain 
that  is  not  evident  during  octamerization. 
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